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Abstract 

We use an Effective Charges (ECH) method to extract A^jg, and hence 
a s (M z ), from event shape moments in e + e~ annihilation. We compare 
these results with ones obtained using standard MS perturbation the- 
ory. The ECH method at NLO is found to perform better than stan- 
dard MS perturbation theory when applied to means of event shape 
observables. For example, when we apply the NLO ECH method to 
(1 - T) we get a s {M z ) = 0.1193 ± 0.0003. However ECH at NNLO is 
found to work less well than ECH at NLO, and the ECH method also 
fails to describe data for higher moments of event shapes. We attempt 
to explain this by considering the ECH /^-function as an asymptotic 
series. We also examine the effect of adding two different models for 
non-perturbative power corrections to the perturbative approximation 
given by the ECH method and MS perturbation theory. Whilst only 
small power corrections are required when using ECH at NLO, it is 
found that these models are insufficient to counteract the undesirable 
behaviour of ECH at NNLO. 

1 Introduction 

In recent years much progress has been made in calculating coefficients for 
higher order terms in QCD perturbative expansions [THS]. This has enabled 
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many QCD observables to be calculated to a greater accuracy [3HS] • These 
more accurate perturbative results can then be used to determine a s {M z ) 
with increased precision [3H8]. However, the result of a calculation performed 
at any fixed order in perturbation theory depends on the renormalisation 
scheme (RS) and renormalisation scale, hr, used. The standard way of deal- 
ing with this dependence is to choose a particular RS (a common choice is 
the modified minimal subtraction scheme, MS [S]) and perform the calcu- 
lation with /ir set equal to some energy, Q. The scale Q is associated with 
the data in some way; for example, when using data from e + e~ collisions it 
is common to choose the centre of mass energy. A theoretical uncertainty 
is then assigned to the result by varying /j,r between Q/2 and 2Q [5|[10|[TTj . 
The choices of \xr and the upper and lower limits of its variation are com- 
pletely arbitrary. This method can also result in theoretical uncertainties 
that are much greater than the experimental uncertainties involved. This 
treatment is therefore somewhat unsatisfactory. 

An alternative approach is to use the Effective Charges (ECH) method 
|12H21j . This method integrates up a /^-function type equation for the ob- 
servable that is analogous to the running equation for the strong coupling 
constant. The resulting integration constant is related to Aqcd ; the charac- 
teristic scale of QCD [5,9.111]. Aqcd 5 known as the dimensional transmuta- 
tion parameter, can be considered to be the truly fundamental parameter of 
QCD [22J. As we shall see, the ECH method directly relates experimental 
observables with this parameter. Usually A^jg- is considered. Once deter- 
mined this underlying universal parameter of QCD can be used to evaluate 
a s (M z ) using the QCD j3- function equation. There is no unphysical depen- 
dence on hr and therefore no need to estimate a theoretical uncertainty to 
account for this dependence. 

In this paper we analyse event shape moments, using both the ECH 
method and standard MS perturbation theory (MS PT), in order to ex- 
tract a s (M z ). The event shape moments have been calculated using data 
taken from e+e" collisions at LEP (DELPHI, OPAL and L3) and PETRA 
(JADE) [23H28] . The DELPHI collaboration performed a similar study [29], 
applying the ECH method and MS PT at next-to-leading order (NLO) to 
event shape means. We replicate this analysis, and also extend the work to 
next-to-NLO (NNLO) using the NNLO perturbative coefficients for event 
shape moments that have recently been calculated by Gehrmann-De Ridder 
et al. in [I]. In a further extension of the DELPHI collaboration's analy- 
sis, we study some higher moments of event shape observables: specifically 
the first three moments of 1 — T, C, Byy, and Y$ as well as the mean 
of pe (see Section [3] for details). As in [2U] we also add power corrections 
to these perturbative models to determine whether the estimation of a s is 
improved. We use the dispersive power correction model, outlined in refer- 
ences [iniEQlET], and also a simple power correction model [29~]l32]. 

The plan of the paper is as follows. We start by outlining the ECH 
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method in Section [2j We then introduce event shape observables and give 
more details of the data used in Section [3j Our results are given in Section 
HI The values of a s (M z ) extracted from the event shape moments using 
the ECH method and standard MS PT, at both NLO and NNLO, are given 
in Section 14,11 In Section 14.21 we show how the ECH method can be used 
to approximate data by working backwards from a fixed value of a s (M z ). 
We provide a commentary of these results in Section 14.31 In Section 14.41 we 
interpret some of the results by looking at the perturbative expansions of 
the event shape moments and values of the A parameter (see also Section [2D . 
Power corrections are then added to the ECH method and MS PT: fits for 
a s (M z ) and a dispersive power correction parameter, ao, are performed and 
the results given in Section [4.51 Fit results for a s {M z ) and a simple power 
correction parameter, kq, are given in Section [4.61 Another commentary is 
given at this point, in Section [4.71 Our conclusions are given in Section 



2 Theory 

Here we shall briefly derive the ECH method |12H21j and show how it can be 
used to extract Ajyjg and hence the strong coupling constant. We shall start 
from the principle of dimensional analysis [22j . Let us consider a dimension- 
less QCD observable which depends on a single energy scale, 1Z(Q). Since 
the observable is dimensionless we have to introduce another, dimensionful, 
quantity A: 



Q, 

where / is some function of the dimensionful quantity A and the physical 
scale Q. We will assume massless quarks; the extension to the case with 
non-zero quark mass has been considered in Refs |17U20j . A must therefore 
be related to Aqcdj the fundamental dimensional transmutation scale of 
QCD. We want to extract A, and to do this we must invert (|2.ip . To find 
the inverse function / _1 we shall start by looking at the derivative of 1Z{Q) 
with respect to Q: 

d1l(Q) = p(K(Q)) (2 2) 

dQ Q { " ' 

We have again used dimensional analysis to determine the form of (|2.2p . It 
can be rewritten as: 

To solve ()2.3p we need to know something about 1Z(Q). We shall take the 
form of its perturbative expansion to be: 

K(Q) = o(l + na + r 2 a 2 + •••), (2.4) 
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where a = (x s ([j,r) / V and p,R is the renormalisation scale. For any observable 
of this type, i.e. one that depends on only a single energy scale, we can 
always ensure that the perturbative expansion takes the form of (|2.4p by 
applying an appropriate scaling and, if necessary, raising to a power. This 
step therefore results in no loss of generality. 

An observable with an expansion of the form of (|2,4p is known as an 
effective charge |2iy33j as it has a running equation like the strong coupling 
does. The renormalisation group equation for a s , rewritten in terms of a, 
is: 

= p {a) = _ 6a 2 (1 + ca + V c n a n ) . (2.5) 

/3(a) denotes the QCD /3-function. The first two coefficients, b= (33— 2Nf)/6 
and c= (153— 19iVy)/12&, are universal. The higher c n are scheme-dependent 

coefficients. In the MS scheme cf 3 = (77139 - 15099N/ + 325A|)/17286 
[SIS]. We can derive the analogous expression for the effective charge by 
setting fiR = Q in (j2.4[) and differentiating with respect to InQ. If we also 
invert (|2.4p to get an expression for a(TZ) and substitute in we find that: 

^| = p{K{Q)) = -bK 2 (l + cK + £ Pn K n ) . (2.6) 

71=2 

Again, b and c are the first two universal coefficients from f|2.5p . The p n are 
specific to the observable TZ(Q); they are RS invariant and Q-independent 
combinations of coefficients from (|2.4p and ()2.5p . For example: 



P2 = T2 + c 2 - ric - . (2.7) 

We now integrate up (|2.6p using asymptotic freedom as a boundary con- 
dition, i.e. requiring that TZ(Q) — > as Q — > 00: 

In — = / —- + K. (2.8) 
Arc Jo P(^) 

This step has introduced an integration constant that we have split into two 
parts: A^. and k. A-ji is a finite dimensionful scale specific to the observable 
1Z. k is an infinite constant which must cancel the singularity in p(x) at 
x = 0. Once this requirement is met there is still some freedom to choose 
the exact form of k. Any variation in the definition of n can be absorbed 
into the definition of Afc. We choose k to be: 

dX (2-9) 



bx 2 (l + cx) 

Adding in this explicit form for k we can now rewrite (|2.8p as: 

Q f°° dx f 71 ^ lb 1 

b In - — = / + / dx 

A n Jn(Q) x 2 (l + cx) J 

= F(K) + G(K) . 



p(x) x 2 (l + cx) 



(2.10) 
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The integration in k has been split into two parts: from to 1Z(Q) and from 
7Z(Q) to oo. Integrating up the first term on the right hand side we get: 



cTZ 



1 + cTZ 



(2.11) 



The form of the second integral, G(1Z), depends on the order at which p(x) is 
truncated. For example, at NLO we have p(x) = —bx 2 (l + cx) which cancels 
with the x Z(i +cx j term to give G(1Z) = 0. At NNLO we have: 



G(K) 



dx- 



P2 



1 + cx)(l + CX + P2X 2 ) 



(2.12) 



This integration has an analytic result with three different forms, depending 
on the relative sizes of the coefficients c and p 2 - 



(i)4p 2 > c 2 ; 
G(ft) = |ln 



A = V4p7 
(1 + cTZ) 2 



l + cK + p 2 Tl 2 



2p 2 



A 



arctan 



2p 2 K + c 



arctan 



(n)4p 2 < c 2 ; 
G(ft) = |ln 



(1 + cKf 



l + cll + p 2 TZ 2 



4p 2 , 
2p 2 



+ 



2A 



In 



2p 2 ^ + c - A 



2p 2 TZ + c + A 



+ In 



c + A 
c - A 



(m) 4/? 2 

G(7e) = 



c 2 ; 



In 



(1 + clZf 



+ (c 2 - 2p 2 



1 



1 



1 + cTZ + P2R 2 
We now rearrange for A-^ in (|2.1U|) : 

An = QT(K)G{K) , 

where 



2p 2 ^ + c c 



F(1Z) = exp 



-1/67?. 



1 



1 W 6 



(2.13) 



(2.14) 



(2.15) 



and G(TZ) = exp (— G(1Z)/b). An is a scheme independent quantity associ- 
ated with the specific effective charge in question. It can be related to the 
universal quantity Aj^r using: 



A 



MS 



(2.16) 



where r = r^ s (pjt = Q). The tilde here indicates that this definition is 
slightly different from the standard one for A^jg (Aqcd in the MS scheme). 
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This is due to our choice of k and, to regain the standard definition of 
, we have to shift n by cln(6/2c) resulting in the relation: 



\is 



2c 



c/b 



2c 



c/b 



-r/b 



(2.17) 



Note that all scheme dependence in A^jg- is contained in the exponential 
factor e~ T l h . Although A^r is a scheme dependent quantity it is possible 
to relate A's in different schemes exactly, providing that the NLO coeffi- 
cients are known in both schemes. One simply notes that, from (|2,16p . 
Aygre r / b = A-r. Since An is scheme invariant we can equate the left-hand 
sides of the expression in two different schemes. For instance for A in the 
MS scheme and A in the MS scheme the relation is: 



„MS 



A 



MS 



exp 



Ams > 



where r MS and r MS are the NLO coefficients in the MS and MS schemes 
respectively. 

From Ajyjg- we can find a s (M z ) using: 



1 



1 + W \ - exp 




7ra 



for NLO 



vr(a + c 2 a 3 ) for NNLO 



(2.18) 



(2.19) 



where W is the Lambert VF-function [36J defined by W(z) exp(W(z)) = z. 
The W-i branch should be used to ensure asymptotic freedom [37] . The cou- 
pling a in (|2.18p is the MS coupling corresponding to a two-loop /3-function 
solution, with the /^-function coefficients c 2 and higher set to zero. This can 
be related to the coupling needed at NNLO, which has a non-zero c 2 but 
vanishing higher coefficients, by an expansion in powers of a. The advantage 
of this technique is that a at two-loops can be written analytically, in terms 
of the Lambert VF-function, and one avoids having to solve transcendental 
equations. 



3 Event Shape Observables 

Event shape observables give information about the topology of an event, 
and are designed to be IR- and UV-safe. They are good observables to 
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use to study a s since they are sensitive to the addition of extra final state 
partons, and hence to extra orders of a s [5]. In this paper we analyse event 
shape moments. The n th moment of an event shape observable, y, is given 

by pannnssj: 

1 f ymax „ da 



(y n ) = — / y n -rd y , (3.1) 

Chad Jo d V 

where y max is the greatest value of the event shape that is kinematically 
allowed and <7had is the hadronic cross-section. 

(y n ) has the following perturbative expansion [HITU]: 



2n l<Av,n+\ 27r J l% n t^oin Q2 



+ (^^) 3 (Cy,n + 2B y , n b In & + A y , n (V In 2 & + 26c In & 
+ . . . , (3.2) 



where b and c are coefficients from the running coupling equation (|2.5p , and 
-4j/,n, By,n and Cy tn are perturbative coefficients that are specific to the event 
shape moment. The NNLO perturbative coefficients (C y ^ n ) for several event 
shapes have recently been calculated [1]. 

In order to apply the ECH method to an event shape moment, we must 
first rearrange (|3.2|) to get it into the form of an effective charge as given 
in (|2.4p . To do this we divide through by A y ^ n and, remembering that 
a = oi s (ijlr)/it, put in appropriate factors of two: 



K{Q) = ^ = a l + ^a + -^ + ... ) . (3.3) 



For clarity we have taken [Xr = Q. Comparing this with (|2.4p we see that 

n = i=- and r 2 = -£^. (3.4) 
2^1 4^1 

In this paper we examine the first three moments of five different event 
shape observables: one minus the thrust (1 — T) [38,39j, the C-parameter 
(C) [USED, the wi de and total jet broadenings (B w and B T ) @21[l3], and 
the three-to-two jet transition parameter in the Durham algorithm (Y3) [44] . 
In addition, we analyse the mean of the heavy jet mass in the E-scheme 
e) 



{p E = {M*/s 

The data come from various LEP experiments including OPAL [23^25] . 
DELPHI [251127] and L3 [28]. Some low energy JADE (PETRA) data are also 
included [24,25]. The data cover a range of energies from 14.0-206.6 GeV. 

All observables are corrected for b and c decays using HERWIG++ sim- 
ulations [MIHT]. Samples of 10 6 events were run with Nf = 3 (only light 
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quarks) and also with Nf = 5 (with b and c quarks). The event shape 
moment was calculated in both of these cases and a ratio of the two quan- 
tities was taken to give a correction factor. This was done for each energy. 
The data in this paper have been multiplied by the appropriate correction 
factors. The correction factors for the means (n = 1 moments) are shown in 
Fig. [TJ Similar trends are seen in the n = 2 and 3 moments, although the 
deviation of the correction factors from 1 tends to get larger at low energies 
for higher moments. These corrections have been compared with those cal- 
culated by the DELPHI collaboration in Fig. 3 of [29J. The factors in Fig. [T] 
show the same energy dependence as in [29] , although our correction factors 
deviate further from 1 at low energies. This is because we have corrected 
for c decays as well as b decays. When only b decays are taken into con- 
sideration the corrections we calculated using HERWIG+-I- agree well with 
those that the DELPHI collaboration calculated using PYTHIA 6.1. 

4 Results 

4.1 Finding A^jg and a s (M z ) from data 

Now that we have corrected data in the form of an effective charge (see 
Section [3]) we can use the ECH method to extract a s (M z ). For every event 
shape moment, and at each energy, a value of A^jg is obtained using (|2. 14|) 
and (|2.16p . An overall value of A^jg is then found by performing a weighted 
average over all the energies for a specific event shape moment. The cor- 
responding values of a s (M z ) are calculated from the averages using f|2. 18|) 
and (I2TT9D . This is done at both NLO and NNLO. 

Table [1] shows the weighted averages of A^jg- obtained from the means 
and the corresponding values of a s (M z ). These results are also plotted in 
the top two panels of Fig. [2j 

For comparison, A^g and a s (M z ) are also extracted from the event shape 
means using MS PT and equation (|3,2p . The results are shown in Table [2] 
and the bottom half of Fig. [2j The first uncertainty on the MS PT results 
is the experimental uncertainty combined in quadrature with the uncer- 
tainties from the perturbative coefficients B y , n and C y ^ n pQ. These small 
uncertainties arise from the numerical integration of the perturbative co- 
efficients associated with the differential cross section ^ in equation (|3.ip . 
The second is a theoretical scale uncertainty obtained by varying [ir between 
fiR = 0.5Q and fj,R = 2Q. The largest difference between the calculations 
using the central value of and the upper and lower variations is taken 
as the uncertainty. The error bars in Fig. [2] represent the combined exper- 
imental and perturbative coefficient uncertainty and the theoretical scale 
uncertainty added in quadrature. In contrast, the uncertainty on the ECH 
results is just the combined experimental and perturbative coefficient un- 
certainty. Results for the n=2 moments are shown in Tables [30 and Fig. [3j 
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Figure 1: Heavy quark mass correction factors for the event shape means 
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and for the n=3 moments in Tables [5]|6] and Fig. [H 



4.2 Approximating data from a fixed value of a s (M z ) 

Another way to implement the ECH method is to work backwards from a 
fixed value of a s (M z ) to get an approximation to the data. In this section we 
fix a s (M z ) at the result extracted from next-to-NNLO (N 3 LO) calculations 
on Z-decays 08]: a s (M z ) = 0.1190 ± 0.0026. This value is close to the 
world average value of a s (M z ) = 0.1184 ± 0.0031 reported in @j. The N 3 LO 
Z-decay value corresponds to = 254 Me V. We can then find An using 
(T2T61) and hence K(Q) by inverting {ZZT3D - At NLO, where G(K) = 1, this 
can be done analytically: 

^NLO(Q) = — f 7 r^7 1 ■ I 4 - 1 ) 



1 + W - exp 



1 + - In - — 

c A^ 

Note that if we rewrite A^ in terms of A^jg- using ()2.16p and compare it 
with (|2.18|) we find that 1Znlo(Q) = &s (mech)A"- The particular choice of 
MS renormalisation scale of ^r, = /Uech = e~ r ^ b Q is equivalent to the ECH 
scheme. This ECH scale corresponds to absorbing radiative corrections into 
the definition of the coupling. 

A NLO ECH approximation to data is then given by (y n )NLO = ^^^^-NLOj 
using (|3.3|) . At NNLO it is not possible to invert (|2.14|) analytically and so 
this is done numerically using Mathematica 7.0 [49]. 

The NLO and NNLO ECH approximations are shown against the data 
in Figs I5TI10I For comparison the MS PT predictions are also shown. These 
are obtained by running the coupling from the fixed a s (M z ) over a range 
of energies using (|2.18p and (|2.19p . at NLO or NNLO as appropriate. The 
data are then approximated using (|3.2p . 



4.3 Commentary on the results 

Since the various tables and figures contain rather a lot of information it 
will be useful at this stage to distill some key conclusions from them. The 
values of a s {M z ) extracted in Section 14.11 agree well with those found by 
the DELPHI collaboration in [29]. They are well grouped, as can be seen 
in Fig. [2] (top left panel); the values for B\y and Y3 are slight outliers. The 
average for NLO ECH agrees very well with the N 3 LO Z-decay calculation 
value [48] . The average for the MS perturbation theory (MS PT) calculation 
at NLO (bottom left panel) sits further from the N 3 LO calculation, and the 
scatter of the central values is larger than for ECH at NLO. The numerical 
values of a s (Mz) are given in Tables [TJ and ^ 

We also demonstrated, in Section 14.21 how the ECH method can be 
inverted to give an approximation to data. Figs. loTfTOl show that the data 
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Observable 


Ams [ M eV] 


a a {M z ) 


Ams [ M eV] 


a s (M z ) 






NLO 




NNLO 


(1-T) 


262 ±4 


0.1193 ±0.0003 


319 ±6 


0.1233 ±0.0004 


(C) 


269 ±4 


0.1198 ±0.0003 


300 ±5 


0.1221 ±0.0003 


(B w ) 


326 ±3 


0.1234 ±0.0002 


312 ±4 


0.1229 ±0.0003 


(Bt) 


265 ±2 


0.1196 ±0.0002 


337 ±5 


0.1244 ±0.0003 


(Y 3 ) 


176 ±5 


0.1126 ±0.0004 


211 ±7 


0.1157 ±0.0006 


(pe) 


268 ±5 


0.1198 ±0.0003 


289 ±7 


0.1214 ±0.0004 



Table 1: Weighted averages for A^jg and the corresponding values of a s (M z ) 
for the means using the ECH method. 



Observable 


Ams [MeV] 


a s (M z ) 


X 2 /dof 


NLO 


(1-T) 


646 ± 9 ± 279 


0.1385±0.0003±0.0107 


52.3/47 


(C) 


597 ± 7 ±247 


0.1365 ±0.0003 ±0.0099 


45.7/41 


(B w ) 


336 ± 3 ±105 


0.1241 ±0.0002 ±0.0062 


39.5/47 


(Bt) 


429 ± 3 ±147 


0.1291 ±0.0001 ±0.0073 


69.7/47 


(Ys) 


268 ± 5 ± 81 


0.1197 ±0.0003 ±0.0056 


61.5/32 


(Pe) 


340 ± 6 ± 90 


0.1243 ±0.0004 ±0.0053 


11.5/14 


NNLO 


(1-T) 


396 ± 7 ± 77 


0.1277±0.0004±0.0041 


51.1/47 


(C) 


378 ± 6 ± 69 


0.1267 ±0.0003 ±0.0038 


44.3/41 


(Bw) 


317±12± 13 


0.1232 ±0.0007 ±0.0008 


37.8/47 


(Bt) 


347 ± 9 ± 36 


0.1250 ±0.0005 ±0.0021 


78.1/47 


(Ys) 


228 ± 5 ± 18 


0.1171 ±0.0004 ±0.0014 


65.1/32 


(Pe) 


295 ± 7 ± 21 


0.1218±0.0005±0.0014 


11.5/14 



Table 2: Aj^g and a s (M z ) for the means using MS PT. The first error is 
experimental and the second is theoretical. 



11 



NLO 



NNLO 



(1 - T) 

(C) ■* - - t 

\Byy) - + * 

(Bt) * * 

(Y 3 ) t» m 

■ ■ ■ ■ i ■ ■ ■ .ii ■ ■ ■ ■ i ■ ■ ■ ■ i ■ ■ ■ ■ I L_i— i—i— Li— i—i— i—i— i—i— i—i 

0.10 0.11 0.12 0.13 0.14 0.15 0.10 0.11 0.12 0.13 0.14 0.15 

a s {M z ) a s (M z ) 



(1-T) - > ► - -~ 

(C) - e 1 - 

(Btf/) - h 1 1 HH 

(Br) - , — i , k+h 

(Is) - - 

I I I I 

0.10 0.11 0.12 0.13 0.14 0.15 0.10 0.11 0.12 0.13 0.14 0.15 

a s {M z ) a s (M z ) 

Figure 2: Scatter plots showing the values of a s {M z ) obtained from the 
means of the event shapes using the ECH method (top two panels) and MS 
PT (bottom two panels). The dotted lines show the value of the coupling 
obtained from N 3 LO calculations on Z decays: a s {M z ) = 0.1190 [3B]. The 
solid lines show the unweighted average for the NLO and NNLO cases. 
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Observable 


Ams [MeV] 


a s (M z ) 


Ams [MeV] 


u s {M z ) 






NLO 




NNLO 


<(1-T) 2 ) 


130 ± 4 


0.1079 ±0.0005 


438 ±22 


0.1298 ±0.0011 


<C 2 > 


123 ± 3 


0.1071 ±0.0004 


466 ±30 


0.1312 ±0.0014 


(B 2 W ) 


166 ± 5 


0.1116 ±0.0005 


203 ±7 


0.1151 ±0.0006 


(B 2 T ) 


0.05 ± 0.02 


0.0520 ±0.0017 


287 ±34 


0.1213±0.0023 


(Yi) 


32 ± 13 


0.0905 ±0.0041 


38 ±16 


0.0924 ±0.0045 



Table 3: A^jg and a s (M z ) for the n=2 moments using the ECH method. 



Observable 


Ams [MeV] 


a s {M z ) 


X 2 /dof 


NLO 


<(1-T) 2 ) 


1047 ±30 ±607 


0.1515±0.0008±0.0172 


23.6/32 


(C 2 ) 


1048 ±29 ±609 


0.1516±0.0008±0.0173 


28.3/26 


(B 2 W ) 


245 ± 9 ± 76 


0.1181 ±0.0007 ±0.0055 


22.5/32 


(B 2 T ) 


775 ±15 ±490 


0.1431 ±0.0005 ±0.0167 


26.9/32 


(Y 3 2 ) 


184 ±29± 63 


0.1132 ±0.0026 ±0.0056 


45.5/17 


NNLO 


((1-T) 2 ) 


467 ±13 ±147 


0.1312±0.0006±0.0070 


23.1/32 


(C 2 ) 


469 ±13 ±147 


0.1313±0.0006±0.0070 


26.4/26 


(B 2 W ) 


207 ± 8 ± 19 


0.1154±0.0006±0.0015 


22.5/32 


(B 2 T ) 


370 ± 7 ±110 


0.1263 ±0.0004 ±0.0061 


26.2/32 


{Y 3 2 ) 


147 ±23± 16 


0.1099 ±0.0024 ±0.0017 


46.9/17 



Table 4: A^r and a s (M z ) for the n=2 moments using MS PT. The first 
error is experimental and the second is theoretical. 
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Figure 3: Scatter plots showing the values of a s {M z ) obtained from the n=2 
moments of the event shapes using the ECH method (top two panels) and 
MS PT (bottom two panels). The solid and dotted lines are as in Fig[2j 
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Observable 


Ams [MeV] 


a s (M z ) 


Ams [MeV] 


u s {M z ) 






NLO 




NNLO 


<(1-T)3) 


90 ± 5 


0.1027 ±0.0007 


647 ±37 


0.1389 ±0.0014 


(C 3 ) 


73 ± 3 


0.0999 ±0.0006 


344 ±57 


0.1248 ±0.0033 


(B 3 W ) 


117 ± 6 


0.1063 ±0.0007 


152 ±7 


0.1104 ±0.0007 


(B 3 T ) 


7.5 ±0.5 


0.0776 ±0.0005 


49 ±8 


0.0953 ±0.0018 


(Yi) 


0.6 ±1.4 


0.0620 ±0.0122 


0.6 ±1.6 


0.0624 ±0.0131 



Table 5: A^g- and a s {M z ) for the n=3 moments using the ECH method 



Observable 


Ams [MeV] 


a s {M z ) 


X 2 /dof 


NLO 


<(1-T)3) 


1249 ± 35 ± 780 


0.1570 ±0.0009 ±0.0200 


7.4/17 


(C 3 ) 


1316 ±33 ±844 


0.1587±0.0008±0.0210 


9.4/17 


(B 3 W ) 


177 ±10± 56 


0.1126±0.0009±0.0051 


8.2/17 


(B 3 T ) 


976 ±21 ±692 


0.1495 ±0.0006 ±0.0205 


9.1/17 


(Y 3 3 ) 


163 ±40± 55 


0.1113±0.0039±0.0053 


67.2/17 


NNLO 


((1-T)») 


493 ±14 ±176 


0.1324±0.0006±0.0081 


7.2/17 




526 ±14 ±191 


0.1340 ±0.0006 ±0.0084 


10.0/17 


<^> 


151 ± 8 ± 13 


0.1103±0.0008±0.0013 


8.5/17 




400 ± 9 ± 141 


0.1279 ±0.0005 ±0.0074 


9.7/17 


(n 3 > 


135 ±33± 13 


0.1086 ±0.0037 ±0.0015 


68.5/17 



Table 6: Aj^g and a s {M z ) for the n=3 moments using MS PT. The first 
error is experimental and the second is theoretical. 
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Figure 4: Scatter plots showing the values of a s (M z ) obtained from the n=3 
moments of the event shapes using the ECH method (top two panels) and 
MS PT (bottom two panels). The solid and dotted lines are as in Fig[2j 
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Figure 5: ECH and MS PT approximations to data at NLO and NNLO for 
((1 - T) n ), with n=l,2,3 from top to bottom. 
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Figure 7: ECH and MS PT approximations to data at NLO and NNLO for 
(Byy), with n=l,2,3 from top to bottom. 
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Figure 8: ECH and MS PT approximations to data at NLO and NNLO for 
{Bj), with n=l,2,3 from top to bottom. 
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Figure 9: ECH and MS PT approximations to data at NLO and NNLO for 
(Y" 3 ra ), with n=l,2,3 from top to bottom. 
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Figure 10: ECH and MS PT approximations to data at NLO and NNLO 
for (p E ). 



for the means are best described by the NLO ECH curves, which lie closer 
to data than both the NLO and NNLO MS PT curves. The exceptions to 
this are (Byy), for which all the perturbative techniques give a very similar 
description of the data, and (Y3), where the MS PT at NLO curve describes 
the data best. 

However the ECH method at NLO works less well for higher moments. 
This is particularly the case for ((1 — T) n ), (C n ) and {BJf) with n = 2, 3, as 
can be seen in the bottom two panels of Figs. [5j [6] and El MS PT also does 
not do a good job at describing the higher moment data. 

ECH at NNLO describes the data less well than ECH at NLO, as shown 
in Section T4.21 The exception to this observation is (I3), in Fig. [9l for which 
the NNLO ECH curve lies closer to the data than that for NLO ECH. Also, 
there is very little difference between any of the curves for (Byy). The 
NNLO ECH curves tend to give a description of the means that is similar to 
or better than the MS PT curves. For the higher moments ECH at NNLO 
is again often very similar to the MS PT curves, although it sometimes fails 
to describe the behaviour of the low energy data as well as MS PT does (e.g. 
(C 2 ) in Fig. E]). 

In the top right panel of Fig. [21 we see that the values of a s (M z ) are 
still fairly well grouped when extracted using ECH at NNLO, with (I3) 
remaining an outlier, but that the average has moved away from the N 3 LO 
result. MS PT at NNLO, shown in the bottom right panel, is seen to have 
improved from the MS PT at NLO case: there are smaller errors due to 
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a reduced scale dependence at higher order, the grouping is improved and 
the average is closer to the N 3 LO calculation. In the higher moment plots 
(Figs. [3] and 2|) the grouping in the ECH method gets worse whereas the 
improvement seen going from NLO to NNLO in MS PT continues to hold. 

These results confirm the success of NLO ECH in providing an accurate 
description of event shape means, as was noted in the analysis of |29j . There 
is generally no need to add a non-perturbative power correction. However, 
for MS PT at NLO additional power corrections are often necessary; these 
are often obtained using a dispersive model of power corrections [TU t lHU l lBT] • 
Whether these power corrections are really required, given the success of 
NLO ECH, is an issue that has been actively debated, in particular at a 
FRIF workshop held in Paris in 2006 [50] . In Section 14.71 we will discuss 
the results of fits using this dispersive model. A new feature of these fits 
is that the ECH method is used in conjunction with the dispersive power 
corrections. 

Before exploring power correction models, we attempt to explain why 
the ECH method at NNLO does not provide a good description of the data, 
even for event shape means. We also try to motivate why NLO ECH does 
not give an accurate description of the higher moments. We do this by 
studying the relative sizes of terms in the asymptotic series expansion of 
p(1Z), and examining the size of An for the variables, as shown in the next 
section. 



4.4 Terms in the expansion of p(JZ) and the size of A^ values 

In order to attempt to explain the results in Section HI we now look at the 
expansions p(TZ). From Section [2] we have: 

dTZ{Q) 



p{K) 



dlnQ 



-bn 2 (l + cR + p 2 TZ 2 + . . . 



(4.2) 



Perturbative expansions in QCD and other quantum field theories are ex- 
pected to be asymptotic series, with an n! growth of coefficients in n th 
order perturbation theory. This divergent large-order behaviour arises from 
renormalon diagrams involving chains of vacuum polarization bubbles (for a 
review see [51]). The large-order behaviour of the p(TZ) function for the in- 
clusive QCD R e + e - ratio has been investigated using renormalon techniques 
in [52], and has been shown to be an asymptotic series. 

A series that is asymptotic to a function, F(a), has the following property 
(see [53] for a discussion of asymptotic series): 



F(a)*J2f n a n if 



n=0 



F(a) 



N 

£ 

n=0 



f n a n 



<\f. 



N+ia 



N+l\ 



(4.3) 



where ~ indicates that the series ]T] / n a n is asymptotic to the function F(a). 
The crucial feature of (|4.3p is that the error associated with terminating the 
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Observable 


n 


cJZ 
(NLO) 


P2K 2 
(NNLO) 


Order at which series 
should be terminated 


((1 - T) n ) 


1 


0.078 


-0.049 


NLO 




2 


0.105 


-0.323 


LO 




3 


0.119 


-0.575 


LO 


(C n ) 


1 


0.076 


-0.030 


NLO 




2 


0.106 


-0.376 


LO 




3 


0.123 


-0.912 


LO 


\ W I 


1 


0.045 


0.004 


NLO 




2 


0.057 


-0.031 


NLO 




3 


0.055 


-0.039 


NLO 




1 


0.066 


-0.060 


NLO 




2 


0.111 


-1.803 


LO 




3 


0.140 


-5.715 


LO 




1 


0.058 


-0.042 


NLO 




2 


0.060 


-0.049 


NLO 




3 


0.060 


-0.053 


NLO 




1 


0.059 


-0.015 


NLO 






Table 7: Size 


of terms in 


-bTl 2 



series at a particular point (after N terms) is bounded by the size of the 
first term neglected (the (N + l) th term in this case). To get the best ap- 
proximation to F(a) we therefore want to terminate the series just before 
the term of smallest absolute magnitude. This is a very important property 
of asymptotic series; although they are divergent, with zero radius of con- 
vergence, one can still obtain accurate approximations by truncating them 
after a suitable number of terms. 

In Table [7] we show the relative sizes of the known terms in f° r 
the event shape moments: clZ, the NLO term, and piR? ■, the NNLO term. 
These numbers are calculated using data from the M z peak. The leading 
order (LO) term in , 1, is universal. The last column indicates the order 
at which the series should be terminated based on the relative magnitudes 
of the currently known terms. Note that, where p^R? is the smallest term, 
it is not possible to say definitively whether the series should be terminated 
at NLO or NNLO without knowing the size of the N 3 LO term. 
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This information can help explain the results discussed in Section 14.31 
For example, the expansions for the higher moments of ((1 — T) n ), C n and 

all suggest that it is best to terminate the series at leading order. This 
indicates that using perturbation theory with this type of expansion, where 
the observable itself is used as the expansion parameter, is inappropriate 
for these event shape moments. As discussed previously, it can be seen 
that applying the ECH method to the higher moments of these observables 
does not work well. However, in some instances (e.g. ((1 — T) 2 ) in Fig. [5]) 
the ECH method at NLO appears to be converging with the data at high 
energies. This implies that we are examining data below the perturbative 
regime of the method for these variables, but that if we go to sufficiently high 
energies the ECH method may be found to describe the data again. As Q 
increases TZ(Q) decreases; for a suitably high energy the NNLO truncation 
will therefore represent the optimal number of terms, and will successfully 
describe the data. Of course this is of theoretical rather than practical 
interest, since we do not have data at indefinitely large energies. 

The expansions for all other variables in Table [7] suggest that it is prefer- 
able to terminate at NLO. As discussed in Section ET3"l it appears that the 
ECH method generally works better at NLO than at NNLO, although a firm 
conclusion cannot be made in these cases since we have no information on 
the size of the N 3 LO term. We can also see that the higher moments of Bw 
and I3 do not display the same undesirable behaviour as higher moments of 
the other variables: see the bottom two panels of Figs. [7] and El 

We also examine the size of A^ for each variable. The values of A-r, 
shown in Tabled are calculated from a fixed value of a s {Mz) = 0.1190 and 
(|2.16p . which when inverted gives = e~ r / b A m . From (gjjty we see that 
7Z(Q) has a Landau pole when Q = A-ji, and will formally diverge there [9]. 
From Table [H the A-ji values for the n=l moments are a few GeV. Therefore 
the data we are using are at energies considerably higher than the Landau 
pole and the ECH method should work well. For the higher moments, the 
values of An increase. We are therefore analysing data in a region of Q 
which contains the Landau pole, and the perturbative ECH method will 
not function until we have data at much larger energies. This is the same 
conclusion that we reached by studying the magnitudes of terms in the p(7Z) 
series. The (Byy) and (Y3 1 ) moments for n = 2,3 are an exception, with 
small values of A^. As stated above, the ECH method performs better for 
these higher moments compared with the other observables. 

4.5 Dispersive power corrections 

In this section we use the ECH method in combination with a dispersive 
model for power corrections [TQHSOJET]. Hadronisation corrections for event 
shape moments are expected to be additive [TUTU]: 

(y n ) = (y n ) P t + (y n )n P , (4.4) 
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n 


((i-m 


(C n ) 










1 


3 


3 


0.1 


2 


1 


1 


2 


14 


15 


1.1 


58 


2 




3 


27 


34 


1.2 


384 


1 





Table 8: Values of K n in GeV 



where (y n ) p t is the perturbative part of the event shape moment and (y n ) n p 
is the non-perturbative part. 

The dispersive model for power corrections accounts for non-perturbative 
behaviour at low energies by replacing the strong coupling constant with an 
effective coupling, a e ff, below an IR cutoff scale, fij. This is done in such a 
way that the integral of the effective coupling up to \xi is finite: 



/ dQa cS (Q 2 ) =oo(mj) • (4-5) 
j o 



1 /w 



This leads to the following non-perturbative contribution to the event shape 
means: 

{y)n P = a y P, (4.6) 

where the a y are numerical factors that depend on the event shape in ques- 
tion |10] . Their values are given in Table [9) Note that since a y 3 = there 
are no dispersive power corrections for (I3™). 
The form of P at NNLO is: 



a,(^) + ^(l + ln^ + !) a2(Mfl) 

+ (46c(l + ln^ + AV 862 (l + ln^ + £) 
fj,i 4bc ) \ hi 2b J 



t2i„ ^ a„ ^ , ^ > \ «s(^) 



+ 46 z In — In 



/-*/ \ W ^ / / 47T 2 

where C F = ^FT ( with ^ = 3 

being the number of colours), and: 
V 18 6 / 9 ' 



Ml 




(4.7) 



245 67 vr 2 11 ll/vr 2 \2\ „ . / 55 



V 108 9 6 3 s / 27 F ' 
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where Ca = N and Np is the number of quark flavours. M is the Milan 
factor, a two-loop enhancement factor [54] . A4 is universal for observables 
with this type of 1 /Q power correction. Given the energy range of the data 
used in this paper, the number of active quark flavours is 5. The form of P 
at NLO (p( JVi °)) is the same as P in (|4.7p except that the a? s term in the 
square brackets is omitted. 

For B\y and Bt there are additional corrections to P. At NLO we get: 



p(NLO) = p(NLO) 
B w 



p(NLO) _ p (NLO) 



TT 



<8C F a s (1 + Ifr) 



3 b_ 

Ka s \ 4 GCp 



TT 



3 b_ 

+ 4 ~ 3C^ 



+ % 



+ 7/0 



(4.8) 
(4.9) 



/40a s (1 + 

where a s (fiR) = a s (e~ 3 ^ i [xr) and t/q = —0.6137. The corrections to P 
for Pvk an d Bt at NNLO have not been explicitly calculated yet but the 
potentially dominant terms are approximated by [10] : 



P 



( 



TT 



Lai 



8^(1 + ^+4, 



+ 



6C F 



(4.10) 



P Bl 



P 



\ 



+ 



Lai 
Irv- 



3 

+ 4 



3C F 



(4.11) 



/ 



In this section we use the forms in (|4.4p and (|4.6p to predict the form of 
the data and then perform simultaneous fits for a s (M z ) and olq. The param- 
eter ao is expected to be universal. For (y) p t we use the ECH prediction at 
NLO and NNLO (see section . In order to facilitate inversion of (|2.14|) 
at NNLO we use Pade Approximant (PA) methods [36,55-58]. We do this 
by writing (gSJ at NNLO as x 2 multiplied by a [1/1] PA: 

/1+ (c— £Z)x \ 
fa) = ~ bx2 1 Twtz + ^(x 3 ) . (4.12) 



We can then integrate up and invert analytically to get: 

^NNLO(Q) = — f 7 T^T ; TTTTxT • ( 4 - 13 ) 



1 



cxp 



1 



In the expressions for the power correction term [ir is set to /Uech = e r l h Q 
(see the discussion below (|4.1|) ) when the ECH perturbative part is used. 
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Event shape 1 - T C B w B T Y 3 p E 

a y 2 3?r 1/2 1 1 



Table 9: The a y coefficients for different event shapes 





Central value 


Upper value 


Lower value 


W [GeV] 


2 


3 


1 


M 


1.49 


1.788 (+20%) 


1.192 (-20%) 


M_r [GeV] 


Q 


2Q 


Q/2 



Table 10: Variation of Ai and 



The IR cutoff should be chosen such that Aqcd <C Hi <C Q. We therefore 
take fij to be 2 GeV. Since the fits for a s and ao have a dependency on the 
choice of cutoff, we vary /i/ between 1 and 3 GeV to assess the extent of 
the dependency. The variation in the fits is used to assign a theoretical 
uncertainty. This is in keeping with the method used in [10]. The Milan 
factor is known at two loop order to be M. = 1.49 ± 20%, so M. is also 
varied around this central value. The central, upper and lower values of 
fij and A4 are shown in Table [101 When the effect of the upward and 
downward variation result in uneven differences the larger of the two is taken 
to contribute towards the theoretical uncertainty. The individual differences 
from the variation of fj,j and M. are added together in quadrature to give 
the theoretical uncertainty. 

Since it was found in the previous sections that generally higher moments 
work less well than the n = 1 moments, we confine ourselves to analysing the 
means in this part of the paper. The results of the dual fits for means at NLO 
and NNLO are shown in Table [TTJ The first uncertainty is the combined 
experimental and perturbative coefficient uncertainty, and the second is the 
theoretical uncertainty (estimated by varying /// and M). These results are 
plotted in Fig. [TTJ We also perform fits using the MS PT prediction at NLO 
and NNLO. These results are shown in Table [T2l and Fig. [T2T In these cases 
Hr is varied (around a central value of /ir = Q), as well as [ii and M, to 
find the theoretical uncertainty (see Table [T0"|) . 

In analogy to Section 14.21 we also perform fits for ao while keeping 
a s (M z ) fixed at 0.1190. This is done for both ECH and standard MS PT at 
NLO and NNLO, and the results are shown in Tables [T3lfT4l and Figs. [T3lfT4T 
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Observable 


a s (M 2 ) 


a 


X 2 /dof 


NLO 


(1-T) 


0.1144 ±0.0007 ±0.0028 


0.436 ±0.011 ±0.006 


52.6/46 


(C) 


0.1133 ±0.0006 ±0.0031 


0.423 ±0.006 ±0.010 


58.5/40 


(%) 


0.1192 ±0.0006 ±0.0014 


0.307 ±0.014 ±0.011 


41.4/46 


(B T ) 


0.1188 ±0.0010 ±0.0038 


0.312 ± 0.027 ±0.053 


59.8/46 




0.1163 ±0.0055 ±0.0010 


0.342 ±0.136 ±0.012 


11.6/13 


NNLO 


(1-T) 


0.1142 ±0.0006 ±0.0017 


0.584 ±0.010 ±0.027 


58.3/46 


(C) 


0.1130 ±0.0005 ±0.0019 


0.539 ±0.008 ±0.026 


67.0/40 


(%) 


0.1184 ±0.0007 ±0.0014 


0.392 ±0.014 ±0.013 


43.1/46 


(B T ) 


0.1184 ±0.0008 ±0.0028 


0.537 ± 0.020 ±0.008 


67.2/46 




0.1168 ±0.0054 ±0.0010 


0.449 ±0.114 ±0.017 


11.6/13 



Table 11: Fits for a s (M z ) and ao using ECH and dispersive power correc- 
tions at NLO and NNLO. 



Observable 


a s (M z ) 


a 


X 2 /dof 


NLO 


(1-T) 


0.1302 ±0.0010 ±0.0088 


0.382 ±0.012 ±0.012 


49.7/46 


(C) 


0.1274 ±0.0008 ±0.0085 


0.364 ±0.007 ±0.013 


48.5/40 


(%) 


0.1191 ±0.0007 ±0.0021 


0.346 ±0.015 ±0.126 


42.7/46 


(B T ) 


0.1252 ±0.0008 ±0.0060 


0.312 ±0.013 ±0.005 


64.3/46 


(pe) 


0.1200 ±0.0060 ±0.0042 


0.323 ±0.138 ±0.018 


11.5/13 


NNLO 


(1-T) 


0.1216 ±0.0009 ±0.0038 


0.427 ±0.011 ±0.020 


49.5/46 


(C) 


0.1197 ±0.0007 ±0.0039 


0.412 ±0.007 ±0.022 


47.9/40 


(B w ) 


0.1179 ±0.0007 ±0.0017 


0.429 ±0.015 ±0.045 


44.3/46 


(Rr) 


0.1212 ±0.0008 ±0.0029 


0.385 ±0.012 ±0.038 


65.1/46 


(pe) 


0.1178 ±0.0056 ±0.0014 


0.389 ±0.127 ±0.033 


11.6/13 



Table 12: Fits for a s (M z ) and «o using MS PT and dispersive power cor- 
rections at NLO and NNLO. 

29 



(1-T) - i ' ' " " m - 

<C> - 1 — - - m 

(Br) - ' ■ - - 

- , , , - - , . 

0.11 0.12 0.2 0.4 0.6 

a s (M z ) a 

1 1 1 1 n 1 1 1 1 1 I 1 1 1 1 1 1 1 1 1 

(1-T) ™ - ^- 

(C) - ~ - ~H - 

(B W ) - l-H 

(Br) - i — h — , - h+H - 

- I H 1 - - I 1 1 - 

0.11 0.12 0.2 0.4 0.6 

a s {M z ) a 

Figure 11: Scatter plots of a s (M z ) and ao from fits of ECH and dispersive 
power corrections for the means at NLO (top two panels) and NNLO (bot- 
tom two panels). The dotted lines on the a s (M z ) plots show the value of 
the coupling obtained from N 3 LO calculations [48J. The solid lines show the 
unweighted averages. 
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Figure 12: Scatter plots of a s (M z ) and ao from fits of MS PT and dispersive 
power corrections at NLO (top two panels) and NNLO (bottom two panels). 
The solid and dotted lines are as in Fig. [TTJ 
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Observable 




a 






NLO 


X 2 /dof 


NNLO 


X 2 /dof 


(1-T) 


0.377 ±0.007 ±0.003 


101.6/47 


0.557 ±0.013 ±0.019 


166.3/47 


(C) 


0.375 ±0.007 ±0.004 


209.6/41 


0.524 ±0.013 ±0.013 


374.8/41 


(B w ) 


0.317 ±0.005 ±0.028 


42.0/47 


0.381 ±0.014 ±0.025 


43.9/47 


(Rr) 


0.314 ±0.004 ±0.007 


59.8/47 


0.526 ±0.017 ±0.036 


68.0/47 


(pe) 


0.282 ±0.010 ±0.008 


11.8/14 


0.403 ±0.015 ±0.016 


11.8/14 



Table 13: Single fits for a , with a s {M z ) = 0.1190 held fixed, for ECH and 
dispersive power corrections at NLO and NNLO. 



Observable 




a 








NLO 


X 2 /dof 


NNLO 


X 2 /dof 


(1-T) 


0.504 ±0.010 ±0.088 


185.3/47 


0.457 ±0.006 ±0.043 


59.8/47 


(C) 


0.433 ±0.007 ±0.063 


198.8/41 


0.418 ±0.005 ±0.028 


49.5/41 


(B w ) 


0.353 ±0.005 ±0.154 


43.0/47 


0.409 ±0.016 ±0.041 


46.9/47 


(Bt) 


0.413 ±0.004 ±0.071 


167.2/47 


0.419 ±0.010 ±0.039 


79.7/47 


(pe) 


0.352 ±0.010 ±0.088 


11.6/14 


0.363 ±0.012 ±0.025 


11.6/14 



Table 14: Single fits for a , with a s (M z ) = 0.1190 held fixed, for MS PT 
and dispersive power corrections at NLO and NNLO. 
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Figure 13: Scatter plots for single ao fits for ECH and dispersive power 
corrections at NLO and NNLO. 
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Figure 14: Scatter plots for single ao fits for MS PT and dispersive power 
corrections at NLO and NNLO. The solid lines are unweighted averages. 
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4.6 Simple power corrections 



For comparison we now use a second model for non-perturbative corrections: 
a simple power correction [29., 32]. For the ECH method this is done by 
adding a Ci/Q power correction (with c\ a constant) to the perturbative 
expansion of 1Z given in (|2.4p . This results in an altered form for the p(TZ) 
function as follows: 



pill) = -bTZ 2 {\ + cR + P2 TZ 2 + ...) + KQ-R-^e- 1 '^ . (4.14) 



This modified form of p{R) is then substituted into (|2.10p . with the pertur- 
bative part truncated at NLO or NNLO, and a fit to data is performed to 
extract A-^. (and hence a s {M z )) and Ko- 

A corresponding power correction can be added to a NLO or NNLO MS 
perturbative part using: 



In this model it is not expected that the fits will yield a universal value of 
Kq. Each particular observable will have a different value of no, unlike the 
dispersive model where the power correction parameter «o is expected to be 
universal. The results are shown in Tables [TMTBI and Figs. [T5TI161 

4.7 Commentary on power correction fits 

In Section T4.5I we examine the effect of adding a dispersive power correction 
to the perturbative model. In Fig. [11] it is seen that the values of a s {M z ) 
and ao extracted from (B\y), (Bt) and (pe) using ECH plus a dispersive 
power correction at NLO agree well with each other. The results for (1 — T) 
and (C) are somewhat discrepant. The values of a s for the other three event 
shapes are not considerably changed from those extracted using pure ECH. 
This implies that only very small power corrections are required for ECH 
at NLO for these observables. At NNLO the extracted values of a s {M z ) 
do not change substantially, but the fitted values of ao become larger. The 
fits for (1 — T) and (C) are still not consistent with those from the other 
observables. The reduced \ 2 values (x 2 divided by the number of degrees 
of freedom) for the fits generally get further from 1 when going from NLO 
to NNLO, indicating that the quality of the fits decreases. 

For NLO MS PT plus a dispersive power correction, shown in Fig. [T2"| 
the error bars are generally larger than for the corresponding ECH fits. This 
is because it is necessary to vary fiR, in addition to the Milan factor and the 
IR cutoff scale, when calculating the theoretical uncertainties. At NNLO 




y,n + • • • 



c 




(4.15) 
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Observable 


a s (M z ) 




X 2 /dof 


NLO 


(1-T) 


0.1230 ±0.0008 


-0.039 ±0.008 


51.3/46 


(C) 


0.1230 ±0.0006 


-0.036 ±0.006 


43.7/40 


(B w ) 


0.1225 ±0.0005 


0.084 ±0.038 


34.2/46 


(&r) 


0.1225 ±0.0004 


-0.054 ±0.008 


60.7/46 


(Ys) 


0.1222 ±0.0021 


-0.237 ±0.045 


31.2/31 


(pe) 


0.1184 ± 0.0052 


0.047 ±0.181 


11.4/13 


NNLO 


(1-T) 


0.1253 ±0.0009 


-0.016 ±0.009 


50.7/46 


(C) 


0.1247 ±0.0008 


-0.022 ±0.008 


43.7/40 


(B w ) 


0.1224 ±0.0006 


0.077 ±0.038 


34.3/46 


(Bt) 


0.1269 ±0.0010 


-0.020 ±0.010 


61.5/46 


(Ys) 


0.1261 ±0.0024 


-0.204 ±0.045 


31.2/31 


(pe) 


0.1196 ± 0.0053 


0.063 ±0.179 


11.4/13 



Table 15: Fits for a s (M z ) and kq using ECH and simple power corrections 
at NLO and NNLO. 



Observable 


a s (M z ) 


K 


xVclof 


NLO 


(1-T) 


0.1374 ±0.0013 ±0.0085 


-0.006 ±0.008 ±0.012 


51.5/46 


(C) 


0.1357 ±0.0009 ±0.0079 


-0.026 ± 0.028 ±0.057 


44.7/40 


{Bw) 


0.1228 ± 0.0005 ±0.0035 


-0.286 ±0.119 ±0.314 


34.6/46 


(Bt) 


0.1305 ±0.0006 ±0.0063 


0.042 ±0.016 ±0.051 


61.5/46 


(Ys) 


0.1260 ±0.0012 ±0.0054 


0.060 ±0.008 ±0.022 


31.6/31 


(pe) 


0.1218 ± 0.0062 ±0.0041 


-0.097 ±0.276 ±0.016 


11.4/13 


NNLO 


(1-T) 


0.1279 ±0.0011 ±0.0025 


0.002 ±0.012 ±0.020 


51.1/46 


(C) 


0.1270 ±0.0008 ±0.0023 


0.018 ± 0.042 ±0.088 


44.1/40 


(B w ) 


0.1223 ±0.0007 ±0.0003 


-0.238 ±0.129 ±0.122 


34.7/46 


(Bt) 


0.1269 ±0.0007 ±0.0013 


0.075 ± 0.020 ±0.033 


61.1/46 


(Ys) 


0.1235 ±0.0012 ±0.0011 


0.074 ±0.009 ±0.010 


31.6/31 


(Pe) 


0.1199 ±0.0060 ±0.0008 


-0.089 ±0.307 ±0.016 


11.4/13 



Table 16: Fits for a s (M z ) and kq using MS PT and simple power corrections 
at NLO and NNLO. 35 
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Figure 15: Scatter plots of a s {M z ) and kq from fits of ECH and simple power 
corrections for the means at NLO (top two panels) and NNLO (bottom two 
panels). The dotted lines on the a s {M z ) plots show the value of the cou- 
pling obtained from N 3 LO calculations [38]. The solid lines are unweighted 
averages. 
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Figure 16: Scatter plots of a s (M z ) and kq from fits of MS PT and sim- 
ple power corrections for the means at NLO (top two panels) and NNLO 
(bottom two panels). The dotted and solid lines are as in Fig. [15] 
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these uncertainties get smaller, as expected, and furthermore the agreement 
between values of a s (M z ) and «o improves. There is no substantial change in 
the x 2 values between the NLO and NNLO fits for MS PT. The description 
of the data by MS PT plus dispersive power corrections is generally better 
than the corresponding ECH fits, as can be seen from the x 2 values in Tables 

unandna 

In Fig. Q21 when a s (M z ) is held fixed and only ao fitted for, good 
agreement is found for ECH at NLO. There is now less discrepancy between 
(1 — T) and (C) and the other three observables. However the x 2 values 
for these two fits are very large compared with the number of degrees of 
freedom, indicating a poor quality of fit. At NNLO the results are more 
scattered and the values of ao required are larger. Again, the reduced x 2 
values for the NNLO fits are generally worse than those for the NLO fits. 

For NLO MS PT, in Fig. [TH there is good agreement in «o but the 
error bars are large. As before the uncertainties decrease at NNLO, as do 
the reduced x 2 values. From the x 2 values in Tables [13] and HH we see that 
at NLO the quality of the fits using the ECH method is comparable with 
or better than those using MS PT. At NNLO MS PT describes the data for 
(1 — T) and (C) much better than the ECH method. 

In Fig. [JS] in Section 14. 6[ when the simple power correction model is 
added to ECH at NLO, very good agreement is found between the extracted 
values of a s (M z ). The fits for kq are generally close to zero, again showing 
that NLO ECH does not require large power corrections; this is an obser- 
vation that was also made in the DELPHI analysis [29]. Slight differences 
in our fit results from the DELPHI collaboration's paper are most proba- 
bly due to differences in the data analysed and the fitting programs used. 
Compared with the results for pure ECH from Fig. [2] we see that adding 
a simple power correction brings the values of a s {M z ) extracted from the 
outliers (Byy) and (Y3) into better agreement, but also increases the un- 
weighted average of a s (Mz)- At NNLO the grouping is less good and the 
average of a s {M z ) is increased even further. 

The fitted values of kq when using MS PT, shown in Fig. [TBI are also 
generally close to zero. The corresponding values of a s {M z ) have therefore 
not altered substantially from when pure MS PT (see Fig. [2]) is used. The 
values of x 2 f° r the ECH method and MS PT are very similar, implying that 
the descriptions of each observable given by the two methods are comparable. 

5 Conclusions 

In this paper we have extended the NLO analysis of e + e~ event shape means 
performed by the DELPHI collaboration in [20], which showed that good 
agreement with experimental data is achieved over a wide range of energies 
by using the ECH approach. We had anticipated that applying the ECH 
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method at NLO to higher moments of event shapes would result in similarly 
good agreement with data, and that extending the analysis to NNLO would 
improve on the NLO results. Surprisingly, it appears that neither of these 
expectations is realised in practice. Higher moments are not described well 
by ECH at NLO, and the NNLO case is even worse. The previously good 
agreement with data for the means deteriorates when the NNLO corrections 
are added to the ECH method. 

In Section 14.41 we attempted to understand the disappointing perfor- 
mance of the ECH method at NNLO by studying the relative magnitudes 
of the successive terms in the ECH /3-function p(TZ) in (|2.6p . as shown in 
Table [7J For the means of event shape moments the NNLO p^R? term was 
the smallest in absolute magnitude, indicating that NLO is the optimum 
order of truncation. Without calculating the next uncalculated p%R? term 
and showing it to be larger in magnitude than the NNLO term, we cannot 
definitively say that NLO is the optimum order of truncation. However, the 
good agreement of the NLO ECH predictions with the experimental data for 
the means supports this interpretation. For the higher moments the NNLO 
term is larger in magnitude than the NLO cJZ term, indicating leading order 
as the optimum order of truncation. The Q evolution of TZ(Q) is controlled 
by the RS invariant dimensionful constant An of (|2.16|) which involves the 
NLO perturbative coefficient r = r^ s (pR = Q). This implies that leading 
order truncation is not physically meaningful. 

Another problem in applying the ECH method to the higher moments 
is the high values of An for these observables. There is a Landau pole 
at Q = An where the integrated /3-function equation diverges. We have 
tabulated the An values in Table For the means the experimental data 
are at values of Q that are much larger than A-ji, which is 0(1) GeV, and so 
applying the ECH method is effective. For higher moments A-ji becomes of 
the same order as the energies analyzed, indicating that the method cannot 
be applied to this data. 

If Q is sufficiently large these problems are avoided as the data will be far 
away from the pole at Q = An. Also, as Q increases TZ(Q) will decrease due 
to asymptotic freedom, and so the optimum order of truncation will become 
NLO or NNLO. It will therefore be possible to apply the ECH analysis 
to the higher moments at sufficiently high energies. This is of conceptual 
rather than practical interest, since at experimentally accessible energies 
the Q-evolution of the higher moments is not adequately described by the 
truncated p(TZ) function. 

As an extension of this work it would be interesting to approximate 
missing higher orders using an all orders renormalon resummation matched 
to the fixed order NNLO result for the p(1Z) function. Such a resummation 
could be done using a dressed gluon exponentiation [59J for instance. 

Although the ECH method at NLO gives good fits to the means at all 
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values of Q we have investigated possible non-perturbative effects using the 
dispersive model of Ref |10p30p31| . which has been widely used in other anal- 
yses of event shapes, and a simple power correction model used in Rcf |29j. 
To include dispersive power corrections in the ECH predictions we have 
taken the renormalisation scale in the non-perturbative power correction 
to be that appropriate to the ECH scheme. When performing fits for NLO 
ECH plus dispersive power corrections we found consistent values of a s (Mz) 
and «o> except for the observables (1 — T) and (C). This is an improvement 
on the fit results gained from NLO MS PT plus power corrections. However 
ECH at NNLO again performs less well than the MS PT counterpart. ECH 
at NLO gave good agreement in a s (Mz) when simple power corrections were 
added. Generally, only small values of kq were required, as expected. We 
conclude that adding power corrections does not improve the performance 
of the ECH method at NNLO, and that non-perturbative power corrections 
are not required by ECH to describe the means at NLO. 

Neither the ECH method or standard MS PT provides a good description 
of the data for higher event shape moments. It is therefore not clear that 
one should be recommended for use over the other. However, the excellent 
results obtained when applying ECH to the means at NLO indicates that 
further study of this area is desirable. 
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